Cosmology with a long-range repulsive force 
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Abstract. We discuss the properties of a cosmology dominated by a charged scalar 
field with a repulsive, long-range self-interaction. The interaction, in the form of a vector 
field with a tiny mass, can have a dramatic effect on the evolution of the universe, with 
interesting consequences - including in some cases accelerated expansion. One charac- 
teristic of the model is an oscillating deceleration parameter, which would potentially 
allow it to be distinguished from other scalar field models such as quintessence. 



1 Introduction 

Recent observations of Type la supernovae with high redshift indicate that the 
universe is entering a phase of cosmological acceleration Identifying the 

causative agent is perhaps the most exciting task for fundamental theory at 
present. There are many candidates. It could be a cosmological constant, the 
need for which was suggested on the basis of other evidence even before the 
supernovae results • Minimally coupled scalars becoming dominant at late 

times was also suggested before the supernovae results f|,|7],|],[§ ■ Since then such 
models have been dubbed "quintessence" [^0) and have received extensive study 
[ pd]Jr^ , [T3| , p^|Jl~^ ] . Nonminimal couplings have also been explored jl6| and recent 
inspiration has been derived from string theory Jl7],[l8| and extra dimensions 
[19|2(|. It has even been suggested that quantum effects may be responsible 



While there is no shortage of models for the physics behind the cosmolog- 
ical expansion, there is no single model which can be considered compelling. 
We have long advocated that there might be interesting cosmological implica- 
tions from long-range forces other than classical gravity ^,^3j. In this paper, 
we investigate the cosmological implications of a universal, repulsive long-range 
force]]. We wish to evaluate the simplest scenario which contains the physics we 

1 For a more detailed discussion and additional references, see Ref. ]24| ]. 
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wish to study: a uniform cosmological scalar with a nonzero local U(l) charge 
density (scalar QED). We find that (perhaps counterintuitively) such universal 
self-repulsion does not in itself create accelerated expansion. However, in combi- 
nation with a suitably chosen potential for the scalar, an accelerating expansion 
with novel observational signatures can be obtained. In particular, the model 
predicts an oscillating deceleration parameter, which could be detected by a 
sufficiently sensitive observation such as SNAPp5||. 



2 Vector long range force 

A simple model for a repulsive cosmological force is a complex scalar field coupled 
to a U(l) gauge field. For example, one might choose a Lagrangian of the form 



<T (D^)* (D u <j>) - \ 9 ap g^F^F pa - V (|0|) 



(1) 



L = C + -m 2 g^A^A„V=9 ■ (2) 



In a homogeneous fluid with a net charge density, the fluid will be self-repulsive. 
However, realizing this situation in an homogeneous, isotropic cosmology is prob- 
lematic. For example, on a closed 3-manifold, the total charge of any infinite- 
range force field must be zero, so the charge density must vanish. (This is a 
consequence of the inability to define the "inside" and "outside" of an arbitrary 
Gaussian surface.) On an open manifold, it is possible to impose a nonzero charge 
density, but only at the expense of isotropy. One must choose a boundary con- 
dition at infinity which selects a direction for the lines of force. These obstacles 
can be evaded by simply making the vector massive pt| , for example with an 
explicit Proca term: 

1 

— 5 
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The current density J M = ie \6 (D ^(p) — (j>* (D M ^>)] is still conserved as a con- 
sequence of global U(l) invariance. The unique solution consistent with homo- 
geneity and isotropy is 

ie {<t)d <j)* - <j)*d Q (f>) 
m 2 + 2e 2 (/>*0 ' [> 

The spatial components of the vector potential vanish. 

The question we wish to ask is: how does a spacetime dominated by such 
a charged scalar evolve? The stress-energy tensor is defined by T^^J—g = 
28S / Sg^ and its nonzero components in this geometry are given by the pressure 
and energy density, 

p = (D^)* (D a <f>) + l -w?Al + V (ct>) , (4) 
p = (D^)* (D a <f>) + l -w?A\ -V (4>) ■ (5) 

We see immediately that the repulsive interaction contributes a term which 
obeys p — p, rather than the p < —p/3 needed for acceleration. This implies 
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that the new term redshifts very rapidly, as a" 6 , and quickly becomes negligible. 
Physically this is because the Universe expands while the mass remains constant, 
so the force eventually becomes short-range on cosmological scales. Therefore, 
the mere presence of a repulsive force does not generically lead to acceleration. 

However, a fully gauge invariant theory will not have a constant mass for the 
vector particle: the mass of the vector will be determined by some scalar field via 
symmetry breaking. The simplest choice is to make the symmetry breaking field 
the charged field <j>, such that m 2 — > 2A 2 </> 2 . The mass of the vector field vanishes 
at <f> — 0, and the force is infinite range. The self-interaction then diverges in 
the massless limit, and we expect interesting cosmological effects if we begin 
with a nonzero charge density and then drive the vector mass towards zero. It 
is convenient to decompose the scalar into a magnitude and a phase, 

0(t) = /(t)e i9 <'>, (6) 

so that A depends only upon the phase, A n cx 9. The energy and pressure 
depend only on the field magnitude /: 



P = f 2 + Jp+V{f) , (7) 
K 



P = f 2 + iff2~V(f) ■ (8) 

Note that the interaction increases the energy density, as one would expect for 
a repulsive force. However, the new term in the energy density and pressure still 
looks like a fluid component with equation of state p — p, and redshifts like 
1/a 6 . One might therefore expect that the interaction becomes negligible at late 
times. This is indeed true when the potential minimum occurs at some nonzero 
value of /. However, the situation is more interesting when the minimum is at 
/ = 0, because then the scalar is prevented from rolling down to its minimum by 
the electromagnetic interaction. For simplicity we will assume a potential with 
a monomial form: 

V(f) = V f b . (9) 
The full set of equations describing the dynamics of the system is 

/ + 3g)/-^ + ^o/- = 0, (10) 



and 



a \ 2 8tt 



3m|j 



K 



f 2 + -^ + Vof 



(11) 



The constant K is related to the charge density J° by: 



(' W) ( "' Vr ' (12) 
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These equations are those of a field moving in an effective potential which de- 
pends upon the scale factor as well as the scalar magnitude, 

V eS (f,a) = V(f) + -^ , (13) 
a J 

and we see that the presence of the interaction gives rise to an "electromagnetic 
barrier" which prevents the scalar field from relaxing to the origin. (This is 
similar to the potential proposed for scalar matter with a global U(l) charge, 
or "spintessence" p7| , p8| .) Due to the dependence upon a the minimum of the 
potential /o is changing as the Universe expands, 

fa oc a^'^ . (14) 

We can use this as a basis for a self-consistent solution to the equations of motion, 
such that the field adiabatically follows the minimum of the potential: 

f(t) oc/ «a«, 

H(t) cx a - 3b /( b + 2 \ (15) 

The solution such that a(t = 0) = is 

o(t) oc ^+ 2 )/ 36 . (16) 

In this case the effective equation of state of the scalar field is determined entirely 
by the exponent b. For example, 6 = 2 corresponds to pressureless dust; b = 4 is 
radiation, p = (l/3)p. However, of particular interest are scalar potentials that 
can provide an accelerated expansion, p < — (l/3)p, which corresponds to b < 1. 

For arbitrary initial conditions, however, the field does not smoothly follow 
the minimum /o of the effective potential. The general solution is a field which 
oscillates about the minimum, with frequency 

W 2 cxa 6 ( 2 ^/( 2+b ). (17) 

The oscillations are driven by the time-dependence of the minimum /o, and do 
not damp with time. Despite the oscillation, however, the behavior ( fi6| ) of the 
scale factor is approximately maintained. Figure 1 shows the evolution of the 
field as a function of scale factor. Figure 2 shows the three components of the 
stress-energy (potential, kinetic, and barrier) as the field evolves. Although the 
potential dominates throughout the evolution, all three terms scale identically 
with a. The precise relationship between the three components is very sensi- 
tive to initial conditions. Note in particular the fact that the kinetic energy 
scales identically with the other terms in the stress-energy, unlike in the case of 
the self-consistent solution /q. Of most interest are, of course, the cosmological 
parameters. Do the small oscillations of the field affect the background? The 
deceleration parameter q follows the behavior of the field (Fig. 3). Note in par- 
ticular that despite the fact that the kinetic energy of the field vanishes when q 
is at its minimum, the equation of state never reaches p = —p due to the con- 
tribution of the barrier term to the stress-energy, and q > —1 at all times. The 
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Fig. 1. Field F = y^8n/3mp l f vs. time for b = 1/2. The dashed line is the solution 
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Fig. 2. The three components of the energy density vs. scale factor for b = 1/2. The 
potential (upper line, blue) and "barrier" (middle line, red) terms dominate the en- 
ergy density, while the kinetic term (lower line, black) is subdominant. All three scale 
identically with a. 
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ln(t) 

Fig. 3. The deceleration parameter q vs. time. The time variable is in the internal units 
of the numerical solution, t = [20, 100] 

system, however is on average in a state of accelerated expansion. Like the decel- 
eration parameter, the Hubble parameter H also exhibits oscillations. Although 
the specifics of the field evolution are model dependent, they arise from the field 
oscillating about the minimum of the effective potential created by the potential 
and barrier terms. So, oscillations can be expected in any model of this type. 
The magnitude of the oscillations on the Hubble diagram can be shown by plot- 
ting the residual in the distance modulus A(m — M) between the exact solution 
and the self-consistent solution, shown in Fig. 4. Figure 5 shows a plot of the 
deceleration parameter q as a function of redshift for the same choice of model 
parameters. The magnitude of the oscillations is A(m — M) ~ 0.04, smaller 
than the accuracy of current supernova la observations, which have errors on 
the order of A(m — M) ~ 0.15 — 0.20(^9). Other choices of initial conditions 
result in a similar order of magnitude for A(m — M), to within a factor of a few. 
Detection of such a signature would present a formidable observational challenge, 
but it is conceivable that future measurements such as those by SNAP[f25[ could 
approach this level of accuracy. We note with interest that there is in fact some 
evidence of small scatter beyond observational uncertainties in existing data, at 
the level of .12 magnitude [p9|. 

3 Summary and Conclusion 

The expansion of the universe is observed to be accelerating. However, what 
physics is responsible for the expansion is at present a completely open question. 
In this paper we look at a cosmological model dominated by a charged scalar 
with a long-range gauge interaction, i.e. a repulsive force. The simplest such 
model, one with an explicitly broken J7(l) gauge interaction, is sufficient to 



Cosmology with a long-range repulsive force 7 



< 




0.01 0.1 1 

z 



Fig. 4. Residual A(m — M) vs. redshift, relative to the self-consistent solution. The 
oscillations are visible at the level of A(m — M) ~ 0.04. 
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Fig. 5. Deceleration parameter q vs. redshift, for reference to Figs. 6 and 7. 

show that a universal repulsive force does not in general lead to accelerated 
expansion. However, the presence of a long-range interaction can substantially 
alter the dynamics of a scalar field. Cosmological acceleration is possible in a 
model in which the mass of the vector field vanishes at the minimum of the 
scalar potential, creating an "electromagnetic barrier" that prevents the field 
from dynamically relaxing to its minimum. Such a model is distinguished by the 
fact that it predicts an oscillating equation of state, which could be observed by 
a sufficiently precise observation of the distance modulus/redshift relation for 
distant supernovae, such as that envisioned for SNAP]25||. 
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